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Abstract 

Due to the long-standing discrepancy between NA35 and NA36 data on A production, 
two drastically different scenarios of strangeness enhancement are still possible. Indepen- 
dent string models, such as the dual parton model, lead to results close to the NA36 data. 
On the contrary, the NA35 results can only be described by introducing full final state 
rescattering of the produced particles. The corresponding predictions for central Pb-Pb 
collisions at CERN energies differ by a factor 3 to 4. Preliminary data on the net proton 
ip-p) rapidity distribution in Pb-Pb collisions favor the independent string scenario. 
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Most of the string model attempts to describe strangeness enhancement in heavy ion 
collisions go beyond the strict framework of independent strings. The modifications consist 
of string fusion [1-8] and, in most cases, final state rescattering of the produced particles 
[3-8]. String fusion is quite natural in string models. In particular, in the dual parton 
model [9] (DPM), string fusion is simulated by qq-qq pairs from the nucleon sea [2,3], 
and, in this way, the string independence is maintained. Final state rescattering, on the 
contrary, is a drastic departure from independent strings models. When fully effective, it 
implies that the results of these models can only be used as an initial condition - the final 
state being strongly modified by the rescattering of the produced particles which drives 
the system towards equilibrium. The model becomes then similar to the so-called fireball 
models - such as hadron gas models and models based on the production of a deconfined 
quark-gluon plasma [10-14]. 

What do we learn from the confrontation of the string models with the available 
data ? A first remark is that in the independent string framework (i.e. without final 
state interaction), the effect of string fusion (or sea qq-qq pairs) is numerically small when 
considering A production. Indeed we know experimentally that the ratio A/ A at mid- 
rapidities is small in heavy ion collisions. In central SS collisions the most accurate value 
of this ratio is 0.24 ± 0.01 [15] and in SW is 0.20 ± 0.01 [16]. We now know that in Pb- 
Pb this value is even smaller : 0.154 ± 0.005 [17]. Since string fusion produces baryons 
in pairs, it is clear that this mechanism can only affect the A yield at 10 % level. The 
situation is quite different in the presence of final state rescattering. In this case A-A pair 
production from string fusion can be much more important provided a sufficient number 
of A's annihilate in collisions with nucleons. 

Based mostly on NA35 results [18], it is widely accepted that a fireball scenario is the 
only possible one - either with QGP formation or with full final state rescattering. The only 
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work where the fireball scenario has been explicitly dismissed [8] , uses, in order to describe 
the data, the VENUS Monte Carlo [7] which actually contains final state rescattering. The 
HIJING Monte Carlo also used by the authors of [8] , which has no final state interaction, 
does not describe either the NA35 [18] or the NA36 data [19] (see fig. lc of ref. [8]). 

In most analysis, the NA36 data on A production are not considered. This is due in 
part to their limited coverage in p±. At present, NA36 data have been fully corrected for 
acceptance, efficiency and decay via unseen channels [19]. (Corrected A and Kq data are 
not given). Moreover, they can be extrapolated to the full p±_ range, or, alternatively, the 
NA35 data can be restricted to the p± range of the NA36 experiment so that a comparison 
of the two data sets is possible. It turns out that there is a discrepancy between them 
which exceeds a factor two [19]. 

The purpose of this letter is to show that the NA36 data on A production can be de- 
scribed in the strict framework of DPM with independent strings. An important ingredient 
in achieving this goal is the novel mechanism of baryon stopping introduced in ref. [20] - 
which increases substantially the A yield at mid-rapidities. No final state rescattering is 
needed. When such rescattering is introduced, the A yield is increased by a factor 2 -h 2.5 
and agreement with the NA35 data is then achieved. It is also shown that the predictions 
of these two scenarios for central Pb Pb collisions are dramatically different and thus a 
clear experimental distinction will be possible. Moreover, recent preliminary data on the 
net proton yield (p-p) in central Pb-Pb collisions [21] give some indirect evidence in favor 
of the independent string scenario. Indeed, the DPM prediction for this difference, with no 
final state interaction [20], is in good agreement with these data. The final state interac- 
tion 7T + N — > K + A, which is mainly responsible for the increase of the A yield, produces 
a corresponding decrease in the N one which destroys the agreement between theory and 
experiment. 

We turn now to the calculation of the A rapidity distribution in pp, central SS and 
central Pb-Pb collisions at CERN energies. We proceed in the framework of the dual 



parton model [3, 9]. In a first step we switch off the qq-qq pairs in the nucleon sea (a = 0). 
Moreover, we also switch off all final state interactions. As discussed in [3] an important 
drawback of the model is a too small baryon stopping. As a consequence, the rapidity 
distributions of both proton and A in central SS collisions have a pronounced dip at mid- 
rapidities not present in the data. This problem has been solved in ref. [20] with the 
introduction of a diquark breaking (DB) component - in which the baryon number follows 
one the valence quarks and is thus produced closer to y* ~ 0. Based on ISR data on p 
and p inclusive production near y* ~ 0, it was concluded [20, 22] that this DB component 
has a cross-section o v £ B = 7 mb, and thus its size is about 20 % - the diquark preserving 
(DP) component, corresponding to the fragmentation of the diquark as a whole, being 
about 80 %. An important result of ref. [20] is that the DB component increases with A 
faster than the ordinary (DP) one. As a consequence, the effect of this mechanism on the 
proton yield in AW and in peripheral SS collisions is rather small - while in central SS 
collisions it grows larger and has the right size needed to reach agreement with the data. 
The prediction for central Pb-Pb is also given in ref. [20]. Recent preliminary data [21] 
nicely confirm both the plateau height and the shape of this prediction. 

With this novel mechanism of baryon stopping the formulae for the A rapidity distri- 
bution in pp and AA collisions are 
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The profile function is given by 

T AA (b) = J d 2 s T A (s) T A (b- s) . 

In numerical calculations, a standard Saxon- Woods form has been used. For central colli- 
sion we take 6 = (taking b < 1 fm the differences are negligeably small) . dN^ L p / dy and 
dNpp L /dy are the ordinary (DP) contributions for leading and non-leading A respectively. 
The corresponding formulae and numerical parameters are given in ref. [3]*. Note that 
with <jdb=o one recovers the expressions of ref. [3]. The DB component is [20] 

= C[p q :(y) + p q :(-y)] , 

n 

P q n v (y) = Z pl(Z)=Z 1 / 2 J z —{l-X-ZT-V 2 . (3) 

This contribution is proportional to the momentum distribution of a valence quark in a 
proton - which in DPM is entirely determined from reggeon intercepts [9, 20]. The constant 
C is determined from the normalization to 0.16. The origin of this value is the following. In 
a NN or pp collision the A average multiplicity is 0.1 [23]. However, in the DB component, 
strangeness production is enhanced by a factor two, since the s-quark can be produced on 
either side of the valence quark carrying the baryon quantum number [20]. Since the size 
of the DB component in pp is about 20 %, we have to decrease by the same amount the A 
normalization in order to keep the average A multiplicity in pp collisions unchanged. This 
leads to the value C — 0.8 x 0.2 = 0.16. Likewise the value of the normalization constant 
a,A in the A fragmentation function given in ref. [3] has to be reduced by 20 % - while the 
corresponding values for A, p and p are unchanged. 

The results for pp — > A are then very similar to the ones obtained in fig. 1 of ref. [3] 
and agree with experiment both in shape and absolute normalization. More precisely, we 
* The non- leading A contribution is proportional to the number of participants, n A , and 
not to the number of collisions, because, at CERN energies, the production of A-A pairs 
in strings which do not involve a diquark is negligeably small. 



aly DB 

dy 



(y) 



6 

obtain < n >a= 0.11 to be compared with the experimental value < n >a= 0.096 ± 0.01 
[23] and (dN/dy) y * =0 = 0.017 to be compared with 0.015 ± 0.005 [24]. 

The numerical results for central SS and central Pb-Pb collisions are given in the first 
columns of Tables 1 and 2. For central SS collisions we see that the result is close to the 
NA36 data. As for central Pb-Pb collisions the prediction for the plateau height is about 
8 A's per unity rapidity. 

We turn next to the effect of the sea qq-qq pairs. The calculation of this effect has 
large numerical uncertainties discussed in ref. [3]. The values given there correspond to 
the maximal possible ones. As discussed above, in the absence of final state interaction, 
the amount of A-A pair production due to this mechanism is limited by the ratio A/A. 
This value is small even at y* ~ [15, 17]. The maximal number of A-A pairs computed in 
ref. [3] would yield much larger values for this ratio. Renormalizing downward this number 
of A-A pairs, in such a way that the experimental ratio A/A is reproduced, we obtain the 
results given in the second columns of Tables 1 and 2. As anticipated, the amount of pair 
production both in SS and PbPb is very small. 

We consider next the fireball or hadron gas scenario by introducing final state interac- 
tion of the produced particles. In ref. [3] it was argued that the reactions mainly responsible 
for the increase and decrease of the A yield are it + N — > K + A and it + A — > K + N, 
respectively. Following [3] the excess of A's is given by 
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(4) 

The particle densities in the rhs of (4) are those obtained in the model without final state 
interaction (their explicit forms are given in [3]) and < a >= 1.5 mb To = 1 fm and 
r = 3 fm. In the calculation one has to divide the in r interval into a large number of 
subintervals (in practice a division into 10 equal subintervals gives a good accuracy). After 
each subinterval one has to evaluate the new values of the p and A densities resulting from 
the final state interaction. These values have to be used as initial conditions for the next 
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subinterval and so on. In order to do so, one has to know the decrease in the proton yield 
associated to the increase, A, in the A one. If the only strange baryon produced by the 
final state interaction were A's, the decrease of the proton yield would be A/2 - with a 
similar decrease in the neutron one. However, since T, ± , ... are also produced (totalizing 
approximately the same excess A) , we assume that the decrease of the proton yield is also 
equal to A [3] . The results for the A rapidity distribution after final state interactions are 
given in columns 4 and 5 of Table 1 for SS and in columns 3 and 4 of Table 2 for Pb-Pb 
central collisions. Again, the first of these columns is the result without sea qq-qq pairs 
(a = in ref. [3]) and the second one is the result with sea qq-qq pairs as computed in ref. 
[3] (a = 0.1)*. 

In central SS collision the A yield has increased by a factor 2 4- 2.5 due to final state 
rescattering and is close to the NA35 data [18]. In central Pb-Pb collisions the final state 
interaction has produced a dramatic increase in the A yield - which is 3 to 4 times larger 
than the corresponding one without final state interaction. In absolute value this difference 
ranges from 15 to 23 units. Such a huge difference should be easy to detect experimentally. 
Data on the A rapidity distribution in central Pb Pb collisions will soon be available. In 
the meantime, it is important to note that there already exists some indirect evidence in 
favor of the independent string scenario (i.e. DPM without final state interaction). As 
already mentioned, preliminary results [21] on the net proton yield (p-p) in central Pb Pb 
collisions are in good agreement with DPM predictions without final state interaction [20] . 
It has been shown above that the latter produces a substantial increase in the A yield and 

* As mentioned above, in the presence of final state interaction, a large A-A pair produc- 
tion as the one obtained in [3] , can be consistent with the small experimental value of the 
ratio A/A due to possible experimental annihilation with nucleons. Since the calculation 
of [3] with (a = 0.1) gives the maximal A-A production which is possible from sea qq-qq 
pairs, the number of A's has to be in between the values given in the two columns a = 
and a = 0.1. 
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a corresponding decrease in the iV one (about 15 units in central Pb Pb at y* ~ 0). The 
decrease in the number of N due ton + N — > K + Ais considerably smaller. Therefore, as a 
consequence of final state rescattering, the p-p yield will decrease destroying the agreement 
between the DPM prediction and experiment. 

In conclusion, it should be stressed that even if forthcoming CERN data confirm the 
DPM prediction without final state interaction (i.e. a A plateau height of about 8 units), 
production of fireballs or QGP droplets remains possible in events rearer than the ones 
considered here. If this is the case it could affect the production of (anti) cascades and 
(anti) omegas. However, the confirmation of the DPM prediction with independent strings 
in central Pb Pb collisions at the level of A production would be quite striking and would 
confine the QGP search to a much lower production level. 
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Table 1 



y* 


(dN/dy) s Jj£ 


(dN/dy) s N s A ^ 




(dN/dy)ff A ^ 





0.72 


0.78 


0.97 ± 0.14 


1.5 


1.9 


2.2 ± 0.3 


0.5 


0.72 


0.77 


0.97 ± 0.12 


1.4 


1.8 


2.1 ± 0.3 


1 


0.71 


0.75 


0.86 ± 0.10 


1.4 


1.7 


2.1 ± 0.3 


1.5 


0.67 


0.70 


0.76 ±0.12* 


1.4 


1.6 


2.2 ± 0.3 


2 


0.61 


0.62 




1.2 


1.3 


1.4 ± 0.2 



A rapidity distribution in central SS collisions at 200 GeV/c per nucleon. Columns 
1 and 2 are the values without final state interactions - respectively without and with sea 
qq-qq pairs (see main text) . Columns 4 and 5 are the corresponding values with final state 
interaction. The NA35 values are read from Fig. lib of [18] and those of NA36 from 
Fig. 14 of [19] with the p± acceptance correction given in eq. (1) (the value with an asterix 
is for y* = 1.25). 

Table 2 



* 
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(dN/dy)2f s b r A 







7.7 


8.4 


23 


31 


0.5 


7.4 


8.1 


22 


30 


l 


6.5 


6.8 


20 


25 


1.5 


4.9 


5.0 


16 


19 


2 


2.9 


2.9 


8.8 


9.1 



Same as Table 1 for central Pb Pb collisions at 160 GeV/c per nucleon. 
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